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Abstract. The deuteron fragmentation by fast protons has been studied both 
near the kinematics of quasi- free proton - proton scattering and far away from 
it. We have concentrated on the interplay between different reaction mecha- 
nisms associated with the antisymmetrization of the initial and final states 
and rescattering contributions. A multiple-scattering-expansion technique has 
been applied to evaluate the reaction amplitude. An essential element of this 
approach in the momentum representation is the use of the effective nucleon- 
nucleon interaction constructed by Love and Franey as a two-body t-matrix 
for the incident proton scattering on a bound nucleoli in the deuteron. Along 
with the five-fold cross sections, the proton analyzing power and the deuteron 
analyzing powers have been calculated as function of the momentum of the 
outgoing fast proton. The results are compared with the data obtained by the 
Gatchina-Saclay collaboration. 

1 Introduction 

Investigations of the proton-deuteron breakup reaction, or more generally of 
nucleon-deuteron reactions, at high collision energies are usually performed with 
the aims of studying the internal deuteron structure, of searching for non- 
nucleonic degrees of freedom, of finding relativistic corrections for the deuteron 
wave function (DWF), etc. To a great extent, they have been stimulated by 
relatively simple recipes of the plane-wave impulse approximation (PWIA) for 
evaluation of the reaction cross sections with unpolarized particles and consider- 
ing beam and target asymmetries with polarized protons and deuterons. One has 
to bear in mind that the approximation in its most popular form is equivalent to 
the so-called spectator reaction mechanism (with the antisymmetrization require- 
ment for three- nucleon states not included). It has already been shown in ref. 
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that the spectator approximation, in which one of the nucleons of the deuteron in- 
teracts with the incoming proton and the other one remains largely undisturbed, 
does not describe the inclusive 1 H(d,p)X data on the reaction cross section and 
the tensor asymmetry T20 for the 0°-angle proton emission at deuteron kinetic 
energies of 1.25 and 2.1 GeV; this is particularly true around the kinematic limit 
where the dominant reaction channel is pd — > ppn. At the same time, according 
to refs. OlSj, the inclusion of more complicated reaction mechanisms, such as 
with the double-scattering (DS) and final-state-interaction (FSI) corrections, to 
the PWIA results has allowed to improve the theoretical description. A similar 
trend has been observed when explaining the angular space asymmetry in the 
deuteron-proton breakup cross sections 4, 5 where one can find an approximate 
estimation of the respective FSI contributions. 

Of course, the understanding of the actual role of two-step reaction mech- 
anisms depends on the reliability of the additional assumptions that are often 
made within the general approach of formal scattering theory, e.g., the Faddeev 
scheme for the construction of transition operators between the initial and fi- 
nal three-nucleon states (see the review [H] and refs. therein). In this context, we 
would like to note that in ref. [B] the author uses just the free nucleon propagator 

go = [E + iO-K}-\ 

with K being the kinetic-energy operator for the three- body system. It enables 
one to perform the calculations without significant complications, due to the 
neglect of the principal value part in the decomposition 

go = P([E-K]- 1 )-i7rS(E-K), 

but it has also to be taken with considerable doubt. 

In the present work we make another endeavour for a proper description 
of the exclusive reaction pd — > ppn that goes beyond the scope of the PWIA 
ansatz whose applicability has no firm justification. Therefore, every time, under 
new kinematic conditions, one has to take care of a coordinated treatment of 
the structure and rescattering effects (in particular, to separate the reaction 
mechanism from possible extraordinary deuteron components). Our approach is 
based on the Faddeev formulation of the multiple-scattering theory for the pd 
system and its full breakup. As in ref. [Hj, we work with the nucleon propagator g$. 
However, unlike ref. [3], we take into account the dependence of the elementary 
nucleon-nucleon amplitudes on the internal motion of a bound nucleon in the 
deuteron, which is affected by the interaction with the projectile proton. 

Along with the Fermi motion influence some off-shell extrapolations are in- 
troduced for these high-energy vertices by using the effective interaction from ref. 
[7j. Special attention is given to their transformation properties under Lorentz 
boosts (in particular, the Wigner rotations). In addition, we avoid seemingly 
artificial truncation in the angular momentum decomposition of high-energy el- 
ementary amplitudes, as made in ref. |3j. 

We apply the matrix inversion method (MIM) 8 , 9 when calculating the half- 
off-energy-shell i-matrix that determines the rescattering in the final (relatively 
slow) p — n pair. The p — n wave functions in the continuum, constructed within 
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the MIM, are convenient from an analytical point of view (cf. the studies of the 
exclusive (e, e'p) reaction on the deuteron in refs. |10Mll| ). Furthermore the spin 
structure of the elementary amplitudes is fully included. This is of particular 
significance for the consideration of various polarization observables. 

The paper is organized as follows. The underlying formalism and the leading 
order terms in the NN t-matrix for obtaining the pd — > ppn reaction amplitude 
are presented in Sec. 2. There one can also find some details concerning the 
application of the MIM. In Sec. 3 we consider the high-energy NN vertices and 
their transformation from the NN center-of-mass system to the deuteron rest 
frame (laboratory system). The polarization observables in question are given in 
Sec. 4. The results of our calculations at two kinematics are discussed in Sec. 5. We 
conclude with Sec. 6. Appendix A contains some useful mathematical relations. 

2 Theoretical Formalism 

The differential cross section of the exclusive reaction 

p(p) + d(0) ^ p( Pl ) + p(p 2 ) + n(p 3 ) (1) 

for unpolarized particles in the laboratory system can be written as 

d 5 a _(2^E p p\ p\ E 2 E 3 \J\ 2 

dpidQ\dfi 2 6 p p 2 E 3 + E 2 (p 2 - pcosG 2 + p\cos0i 2 y 

where 2 is the emission angle for the outgoing slow proton and 0\ 2 is the angle 
between the momenta pi and p 2 of the fast and the slow protons in the final 
state (see Fig.l). 

We introduce the orthonormal basis 



P PXPl /q\ 

z = — y = r , x = y x z (3) 

\p\ \pxpi\ 

assuming that p and p\ determine the reaction plane so that y is perpendicular 
to it. 

Following the Alt-Grassberger-Sandhas (AGS) formalism of three-body col- 
lision theory (see, for instance, refs. ^21E])> the reaction amplitude J is given 
by the matrix element 

W = ^ < 123 lt 1 " (M) " (l,3)]Z7 i|l(23) >= 5(p - Pl - p 2 - p 3 )J (4) 

taken on the energy shell: 

E = E p + m d = E 1 + E 2 + E 3 . (5) 

Here E is the total energy of the three-nucleon system, the deuteron mass, and 
(i, j) the permutation operator for two nucleons. The state 1 1 (23) > corresponds 
to the configuration of nucleons 2 and 3 forming the deuteron state and nucleon 
1 being the projectile, whereas the state < 123| represents the free motion of 
three nucleons after the reaction and is the product of single nucleon states. 
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Iterating the AGS equations for the rearrangement operators, one obtains 
the transition operator Uqi as 



01 



1 + *2 + h + hg t 2 + hg t 3 + t 2 g t 3 + hg Q t 2 + 0(t 3 



(6) 



By definition, t 3 = t\ 2 is a two-nucleon transition operator (the other can be 
obtained via cyclic permutations). Obviously, the term g^ 1 does not contribute 
to the on-shell amplitude (|IJ). 

We are interested in such kinematic situations where, firstly, one of the final 
protons carries away a considerable portion of the initial momentum p and, sec- 
ondly, the momentum q = \p—pi\, transferred to the p-n pair, is small in compar- 
ison with the nucleon mass m. It enables us to use various nonrelativistic models 
for the DWF. Under such conditions we retain the first- and second- order contri- 
butions to the multiple-scattering expansion (JBJ), neglecting the corrections due 
to the so-called recoil reaction mechanism for which the high-momentum proton 
leaves the deuteron without direct knock-out. In addition, small contributions 
from the DS with participation of a fast nucleon, i.e. the t(fast)got(fast) terms, 
are disregarded, what is in accordance with the conclusions in ref. After this, 
the matrix element Upd-* ppn can be approximated by 



U p d- 



•ppn 



v / 2< 123|[1 - (2,3)][1 +i 23 5o]^ m |l(23) > 



(7) 



with the symmetrized NN t-operator tV^ 1 = [1 — (l,2)]ti2- 

The three-body operator uj 23 = 1 + £2350 acting on the state < 123| yields 



< 123|w23 = hni <123|- — - 1 + t 23 

e^0+ E + l£ - K 

= lim < 123|- — ^ — , 



1 



E + ie-K 



(8) 



where V 2 % is the potential operator of the two-nucleon interaction . Here we 
have employed the relation 



1 



1 



E + ie-K -V^ 



+ 



1 



*23" 



1 



2:! 



E + ie-K E + ie-K E + ie-K 



(9) 



for the three-body operator t 2 %, which obeys a Lippmann-Schwinger (LS) equa- 
tion with V23 as driving term. 

Further, since K = K\ + K 2 + K 3 and, by convention, i>Q|123 >= £^123 >, 
one can rewrite Eq. Q in the form 



< 123|cj 23 



lim < 1231 — 

e^0+ E 



IE 



Et + ie-K 2 -K 3 - V 2 



23 



(10) 



or 



< 123|w 23 =< 123| [1 + t 23 (E - E^g^E - EJ] 



A possible way for the construction of a NN interaction starting from a fundamental Hamil- 
tonian of interacting meson and nucleon fields has been shown recently in ref. |14| 
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The operator g 23 {E — E\) = (E — E\ + iO — K% — K^)' 1 is a free propagator 
for the (23)-subsystem, and the scattering operator t2 3 (E — E\) satisfies the LS 
equation 



t2z{E — E\) — V23 + V23523 {E — E\)t2z{E — E\). 



(12) 



Let us rewrite the matrix element (0 indicating explicitly the particle quan- 
tum numbers, 



U, 



pd-^ppn 



V2 < Pl/XlTi,p 2 A t 2T2,P3^37"3|[l ~ (2, 3)]u 23 |PM T ) V'lAfrfOO (23) >, 



with the spin and isospin projections denoted as fi and r respectively. We assume 
that t = ti = T2 = 1/2 and T3 = —1/2. Inserting the unity operator 



dp'\p' h't' X p'fjfr'], 
we arrive at the following expression for the reaction amplitude 
J = ^C TT , < -ti 2 -fj, 3 \SMs >< -ti 2 -»' 3 \SM' s > 

t ST (E — Ei) 



TT 



dp r '(p r ,SM s 



1 + m 



Pr — Pr + ^0 



pAsm's 



< P1H1, (Pr + 9/2) v' 2 \t% m (E - E' 3 )\pn, {p' r - q/2) //' > 

< fJ,"^iM d {pr' ~ 9/2) > "(2^3), 



(13) 



where E' 3 = y m 2 N + (p/ — q/2) 2 . Here we have introduced the momentum 

transfer q = p — Pi, relative momenta p r = |(p2 — P3) and p' r = \{p' 2 — P3), 
and isotopic coefficients 



Ctt 1 = 2 



-5t'o + 1 + 



-I) 5 



(14) 



Henceforth, all summations over dummy discrete indices are implied. 

In momentum representation the DWF i^iM d (k) with spin projection is 
written as 



|VW*0>=E E <LM L lM s \\M d >u L (k)Y^(k)\lM s >, (15) 



L=0,2 M L =-L 



with the spherical harmonics Yr L (k) and the Clebsh-Gordon coefficions in stan- 
dard form. In our calculations, we have employed the following parametrizations 
of the S- and D- state wave functions 



uo(p) 



c, 



tt i-r 1 a 2 + p 2 



u 2 (p) 



2 >p Dj 



(16) 



as proposed in refs. ^3E3 • 
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The wave function of the final p-n pair 

/,(-) 

J p r SMsTM T 



< ^Jsm,tmM'SM' s TM t >= S( Pr -p' r )S Ms M' + 



-.2 _ n l2 



contains the FSI part, which can be taken in different ways. First of all, we note 
that in ref. [I] the t-matrix elements were obtained by solving Eq. (|12|) with a 
separable NN interaction. Of course, these solutions are not able to reproduce 
the off-shell behaviour of the NN t-matrix with more realistic two-body forces. 
Another approximation relies on the neglect of the principal value part in 



-mb(vl ~ P?) + P 



1 



Then, the FSI term can be expressed through the on-energy-shell t-matrix, i.e. 
the respective phase shifts. Third procedure ( see, ref. ^7] and refs. therein) 
uses some artificial off-shell extrapolation of the t-matrix with energy-dependent 
prescriptions for the introduction of a proper form factor. 

In this paper we use the MIM applied to study the deuteron electrodisinte- 
gration jlUUllj . As in ref. ^U], we consider the truncated partial- wave expansion, 

< V£;L S TM T \Pr'SM' s TM T >= 5 MsM!s 5( Pr - Pr ') + (19) 

J max J 

£ E ^Pr) < lmSM s \JMj > V^(Pr') < l'm'SM' s \JMj > Y^'iPr'), 

J=0 Mj=-J 

where J max is the maximum value of the total angular momentum in n-p partial 
waves and a = {J,S,T} is the set of conserved quantum numbers. The radial 
functions ij)f v (p/) are related via 

HFiPr) = £ OuwW) ~ 6{P ' r ~ Pr) Sll> (20) 
yf Pr 

to the partial wave functions <£>p Z ,(p[.), which have the asymptotics of standing 
waves. The coefficients Oun can be expressed in terms of the corresponding phase 
shifts and mixing parameters jlUj . 

Within the MIM, the functions ^{p^) can be represented as 



where the coefficients Bfi,(j) fulfill a set of linear algebraic equations approxi- 
mately equivalent to the LS integral equation for the n — p scattering problem. 
Here N is the dimension of this set, pj are the grid points associated with the 
Gaussian nodes over the interval [-1,1] and pn+i = Po (details can be found in 
ref. |18j). It should be noted that in this way the nucleon wave function is ex- 
pressed by a series of 5-functions allowing one to reduce a triple integral in Eq. 
(fT3|) to a double one. In addition, the method offers the opportunity to consider 
the nucleon wave function in the continuum directly in momentum space what 
simplifies all the calculations a lot. 
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3 High-Energy Vertex and Its Transformations 

As seen from Eq. (13), we have to deal with the high-energy matrix elements 

t T (E on - s heir,Pl^l,P2^P^l:P'^2) = 
< PlVl,P2tf2\ tT ' { E on-shell)\Pll>\,P 7*2 >> ( 22 ) 

where p\ ~ p > m and E on _ s f ie n > 2m , which are determined not only on the 
energy shell 

E Pl + E p > 2 = E on _ she u = E p + Epi (23) 

but largely beyond it for N-N scattering with arbitrary values of the total mo- 
mentum p + p ' = pi + P2 ' (see Fig. 1). The latter is conserved owing to the 
requirement of translational invariance. 

The frame dependence of these off-shell quantities, in particular, their expres- 
sion in terms of the t-matrix elements in the N-N center-of-mass system (c.m.s.) 
can be found making use of some prescriptions |19l I2()j of the relativistic po- 
tential theory (RPT) for relating a two-body off-energy-shell f-matrix between 
different frames of reference. First, let us recall the equation 

tiE-k'^k'^M^iM^) =N'N < /iil^V^'i)!^ > (24) 

< n' 2 \W\v! ,k' 2 )\v' 2 > t(E,K-k'*v[v' 2 -k*vi,V2) 

< u^W^k^lm >< u 2 \W(u,k 2 )\fJ,2 >, 

which connects the on-momentum-shell matrix elements in two bases 2 . One of 
them, {\k\fJLi, k 2 fi 2 >}, is a direct product > ®\k 2 fi 2 >} where the single- 

nucleon states \k a fj, a > satisfy the following orthonormality and completeness 
conditions 

< K^'a\ k aHa >= 6(k' a - k a )5^ a (25) 

and 

J dk a \k a fi a >< k a n a \ = 1. (26) 

By assumption, they are related to the states |[m, // >= \^2Ei c \k^i > of 
the so-called canonical basis, which is well known in elementary particle physics 
(see, e.g., Ch. IV of ref. |21|). The latter is transformed by an irreducible repre- 
sentation of the Lorentz group. In particular, we have under the Lorentz boost 
L(v) corresponding to the motion of the reference frame with velocity v, 

1 1 

\[m,-],k a ,ij a > tr =\[m,-},k t a r ,u a ><u a \W(u,k a )\n a > (a = 1,2), (27) 

where 

W(u, k a ) = exp(-kr a • u a ) = 

/.. . °a ■ [U X fc a ] \ 

cos — si — i-7z ; 7^ — > (28) 

2_i (1 + u )m + Ef + E a J 1 ; 

2 Henceforth, the asterisk is used to mark quantities in the c.m.s. 
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is the Wigner rotation operator in the spin space of the a-th nucleon. Here we 
introduce the four-velocity u = (uq,u) = (uo,uqv) with a Lorentz scalar u 2 = 1 
and denote the nucleon energies in the initial and moving frames of reference by 
E a = \J k 2 + m 2 and E l J = \J k^ 2 + m 2 . The angle uj a of this rotation about 
the axis [u x k a ] is determined by the relation 

tan^ = 7 l[UXk H (29) 

(see, for instance, Eq. (32) in ref. [22] or Eq. (C.7) in ref. ). In formulae (|27j) - 
(|29|) the symbol tr belongs to the transformed quantities, so that the transformed 
four-momentum is 

kl r = (Ei r , kt) = L(v)k a = L(v)(E a , k a ). (30) 

Since we are interested in a specific transition to the N-N c.m.s. simulta- 
neously for the initial and final N-N states, all we need is to put in Eq. (|24jl 
u = ~K I ^[s and u' = K' j\fs' \ where K = k\ + &2 = k[ + k' 2 = K' is the total 
momentum of the initial(final) N-N pair. Normally, the Lorentz scalars s and s' 
are defined as s = (h + k 2 ) 2 = E* 2 and s' = (k[ + k' 2 ) 2 = E*' 2 . 

The corresponding basis {\Kk* nifi2 >} is composed of two-nucleon states 
with relative momentum k* between two nucleons, measured in their c.m.s., and 
total momentum K. In accordance with the customary recipes (see, for instance, 
Sect. 28 in ref. US])) every state \Kk*fiifi2 > can be represented as a product 
state | [m, |], > | [to, |], —k*,[i2 > which is subject to a subsequent Lorentz 
boost with velocity u = —K/E*. Doing so, we get 

\Kk*mn2 >= N\k lVl > \k 2 u 2 >< v^W^K/E^kx^m > 

<u 2 \W\K/E*,k 2 )\fi2> ■ (31) 

We have used the property W(—u,k) = W'(u, k) and the relations 
W(K/E*,k*) = W(K/E*,ki) and W(K/E*,-k*) = W(K/E*,k 2 ) . In this 
context, note that for given momenta k± and k2 in an arbitrary frame of reference 

k* = kl = -k* = E2 + E * 2 k, - El+El fc 2 , (32) 
1 2 E 1 + E 2 + E* E 1 + E 2 + E* y ' 

where E* = E\ + E 2 is the respective invariant mass (the total energy in the 
c.m.s.). 

In order to provide the orthonormality and completeness relations 



< K'k'y'^Kk*^ >= S(K' - K)S(k'* - k*)8^^ >2 (:$:$) 

and 



dKdk*\Kk* >< Kk*twjL 2 \ = 1 (34) 
the normalization factor N in Eq. (|31|) (see also Eq. (|24jl ) should be 

iV=[J(fei,fe 2 )]-5 ; (35) 
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where J(ki,k 2 ) is the Jacobian for the transformation between the kik 2 - and 
Kk* - representations, i.e. 

1(1 1 \ ^1+^2 ^1-^2 (tl „s 

J(k 1 ,k 2 )= 1 (36) 



Analogously, one has 



E'* -\- E f * E f E f 
J(k 1 ,k 2 )= ( 37 ) 



and the normalization factor N' for the final N-N pair is 

N' =[,J(k[,k' 2 )]-K (38) 

The next result [20] 



t(E; K, k'*u[u' 2 - k*ui, v 2 ) = Ft(\/ E 2 - K 2 ; k'*v[v f 2 ; k*v u v 2 ), (39) 

with 

E VE 2 — K 2 + y/(E[ + E> 2 ) 2 - K 2 



y/E 2 - K 2 E + E[ + E' 2 

^E 2 - K 2 + y/(Ei + E 2 ) 2 - K 2 
E + Ei + -E'2 



(40) 



enables one to express the fully-off-shell matrix elements with nonzero total 
momentum through the others with zero total momentum. The latter are the 
fully-off-shell i-matrix elements in the c.m.s.. 

Using Eqs. (|24[) and ()39|) in our evaluation of the matrix elements ()22|) . we 
approximate them with the ones in the N-N laboratory system (LAB) with p ' = 
0. This appears to be reasonable since p ' = p' r — q/2 (see Eq. (|13|0 is an internal 
nucleon momentum in the deuteron. In addition, we neglect the deuteron binding 
energy, and thus 



Eon-shell — E p + m = E, E' — E[ + E 2 — E Pl + Ep> 2 , 

s = 2mE, s' = E' 2 -p 2 . 

Of course, in order to avoid any confusion, one should distinguish between the 
total two-body energy E in this section and the collision energy in Eq. ©. 

Then the transformation of interest for this half-off-shell t-matrix can be 
written in the compact form, 

t LAB ^t T \E-pi^q^p^^ 2 )=<PW\u\pi)W\u',q)t CM (41) 



with 
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or 



T P \ fi T pi \ I EE ' + ^ 



where T p (T Pl ) is the kinetic energy of the incident (outgoing) proton, and t CM = 
k'* ; k*) is the half-off-shell i-matrix in the N-N c.m.s. Further, a simple 
calculation yields 



W t (w',pi) = cosy[l + io-i- ? /tany] ! (44) 



W\u',q) = cosy \\-%a 2 -y tan (45) 



tan Y = 7pK- 1 smeu (46) 



tan T = — — sm^, (47) 



and 



with 



and 



where 6\ is the angle between the momenta p and p\, i.e. the emission angle of 
the fast proton in the deuteron rest frame. Explicit expressions for the coefficients 
K\ and K 2 are given in the Appendix. From Eqs. (|44|) and (|45|) we see that the 
Wigner operators correspond to rotations in opposite directions about the y - 
axis normal to the reaction plane. 

When evaluating the half-off-shell elements t{yfs] k'*; k*) in Eq. Q41j) we have 
started with the Love-Franey model used to construct the on-shell N-N t-matrix 
in each N-N channel (triplet-odd, triplet-even, etc.). According to ref. [7j (see 
Eq. (12) in the first paper), the corresponding matrix elements are expressed 
through the effective N-N interaction operators sandwiched between the initial 
and final plane- wave states. We extend this construction to the off-shell case using 
the same operators with complex strengths depending on ^/s and allowing the 
initial and final momenta to assume the current values of k* and k'* . Obviously, 
such an off-shell extrapolation does not change the general spin structure 

t%% = A + B(cr 1 n*)(cr 2 n*)+C(cr 1 + cr 2 ).n* + 
DfaDfaf) + F(<7iQ*)(<t 2 Q*), 



where 



Q* = , n* = (48) 



inherent in the on-shell N-N i-matrix in c.m.s. (see, for instance, ref. |24j). Of 
course, this structure is subject to certain changes under the transformation (|41|) 
so that 

t% b N = cos y cos y#{F x + F ^iy) + F 3 (a 2 y) + F 4 (a lX )(<T 2 x) + 

F 5 (a iy )(a 2 y) + F & (<T lZ )(cr 2 z) + F 7 {a lX )(a 2 z) + F%(<t iZ ){<t 2 x)} . (49) 

Relationships between the structure functions Fi in the laboratory and A to F 
in the center- of -mass systems are given in the Appendix. 
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4 Polarization Observables for Proton-Deuteron Breakup. 

The reaction amplitude J resulting from our calculation, allows us to find the 
spin polarization coefficients 

Ci,j,k,l,n = -^rTr {jViajJ^G k Gi(T n } (50) 

with 

C = Tr{jJ^, (51) 

where the first and second indices i and j refer to the deuteron and the incident 
proton polarizations, respectively, and the remaining ones specify the nucleon 
polarizations in the final state. The standard operators = 0,1,..., 9) intro- 
duced by Goldfarb (see, e.g., Eqs. (2.17) in ref. j2Sj) form the basis for the 3x3 
spin 1 space. If a particle in the entrance channel is unpolarized or the polariza- 
tion of a final nucleon is not analyzed, the corresponding index is assumed equal 
to zero. 

In this paper, we consider the first- order observables related to the polariza- 
tion of initial particles for coplanar geometry. In this context, let us determine 
the analyzing power of the projectile proton 

Co l2 /,o,o,o = 7T ImS )Z < ^1,^2, m 3 \J\-M d >< M d , --\J\mi, m 2 , m 3 >*, (52) 
Oo — £ I 

and the vector and tensor analyzing powers of the deuteron 
a/2 

^,0,0,0,0 = —ImS^{< m 1 ,m 2 ,m 3 \J\ml >< 0m\J\m 1 ,m 2 ,m 3 >* - 
^o 

< mi, m 2 , m 3 \J\m — 1 >< 0m\J\mi, m 2 , m 3 >*}, 
C w ,o,o,o,o = yZi\ < m i> rn 2 ,m 3 \J\ml > | 2 + | < m 1 ,m2,m 3 \J\m - 1 > 

2| < mi,m2,m 3 \J\mO > | 2 + 

6Re < nii,m2,m 3 \J\ml >< — lm\J\mi, ni2, m 3 >*} 
Czz,o, 0,0.0 = VV| < m 1 ,m2,m 3 \J\ml > | 2 + | < m 1 ,m2,m 3 \J\m - 1 > | 2 
2| < mi,m2,m 3 \J\mO > | 2 }. 

Of course, the symbol ^ means summation over the particle spin projections. 
Here, unlike in Eq. ()13|) (more exactly, in its l.h.s.), these projections are explicitly 
indicated. As a consequence of parity conservation, the vector analyzing powers 
of the deuteron C x>0)0l o,o, C 2i0 ,o,o,o and the proton Co^q.o.o, C , 2 ,o,o,o are equal 
to zero in coplanar geometry. Moreover, the tensor analyzing power C XX! o, o,o,o 
is not independent owing to the relation C XX! o,o,o,o + Cro,o,o,o,o + C 22j o,o,o,o = 
(cf. Eq. (2.18) of ref. [2Hj). Thus, we have only the above four independent first- 
order observables. 
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5 Results and Discussion 

Using expression (fT3)) for the reaction amplitude, the exclusive cross sections 
and the first-order polarization observables have been calculated at a projectile 
energy of 1 GeV in the deuteron rest frame. All the calculations have been per- 
formed with the Paris potential JHj • We have confined ourselves to the coplanar 
geometry (with momentum p2 lying in the reaction plane) and considered two 
kinematics. 

Kinematics I is equivalent to the Gatchina experiment |17| I2(ij as the pp 
coincidence cross sections have been measured on the left of the quasi-free peak 
(QFP), far from its maximum (Fig. 2). Recall that the latter corresponds to 
proton elastic scattering on a free proton at rest in case of a spectator momentum 
P3 = 0. In fact, the four-momentum conservation for reaction (^Q) in the deuteron 
rest frame yields the Bjorken variable x = q^/i^lmuj) = 1, where q 2 = q 2 — u> 2 = 
—t. Of course, the equality x = 1 is valid only up to small corrections due 
to binding effects in the deuteron. The FSI distortions have been taken into 
account in p — n partial scattering states with J max =2,3, and 4. We see that the 
results obtained for the case with J m ax = 3 are practically converged and may 
be assigned a reasonable predictive power (cf. Figs. 2-6). 

The experimental conditions have been chosen so that FSI contributions 
could be neglected. Indeed, the FSI give insignificant corrections for the dif- 
ferential cross section with respect to the PWIA predictions (Fig. 2). However, 
their effects are more pronounced for the polarization observables. In this con- 
nection, especially the tensor analyzing power of the deuteron C TO ,o,o,o,o (Fig. 
5) and the vector analyzing power of the proton 6*0^,0,0,0 (Fig. 3) are of great 
interest. Note in particular that the PWIA curve for C yy ,o,o,o,o (Fig.5) has a very 
quickly changing behaviour and comes from a negative region up to nearly the 
upper limit of 1 while the full result lies only in the positive region, has a modest 
slope and varies from ~ 0.5 to 0.8. We observe a different behaviour for the other 
tensor analyzing power of the deuteron C 22 ,o,o,o,o (Fig. 6). One can see that this 
observable in PWIA assumes small negative values while the calculations taking 
account of the FSI produce a positive result with very quickly changing values. In 
addition, the result obtained with J max = 2 is significantly different from the one 
with J m ax = 3, what is not so evident for the other observables. For the vector 
analyzing powers of the proton 60,^,0,0,0 (Fig. 3) and the deuteron C y ,o,o,o,o (Fig. 
4) the PWIA and FSI curves have a similar behaviour, but their absolute values 
are about two to three times different. Thus, we can conclude that the polariza- 
tion observables are very sensitive to reaction mechanisms. Unfortunately, only 
the differential cross section has been measured under such kinematic conditions. 

The same observables in kinematics II are presented in Figs. 7-11. The pro- 
jectile proton energy is 1 GeV as in the previous case. The scattering angles of 
the fast and slow protons are 6\ = 6.77° and 62 = 83°, respectively. In such 
kinematics we have more favourable conditions for the applicability of our ap- 
proach. Really, while the q (p^) - values for kinematics I range from 500 (300) 
to 650 (570) MeV/c, the momentum transfer q and the neutron momentum p^ 
for the kinematics II with a cut at p\ = 1.60 GeV/c do not exceed, respectively, 
Qmax = 217 MeV/c and p™ ax = 180 MeV/c which is in complete accordance with 
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the underlying assumptions. 

The FSI corrections calculated for J max =1 and 2, are denoted by the dash- 
dotted and solid lines, respectively. The agreement between them is quite sat- 
isfactory, especially for the differential cross section (Fig. 7), where the lines 
for Jmax =1 and 2 are indistinguishable. As before, the PWIA predictions are 
denoted by dashed lines in all the figures (Figs. 7- 11). 

Unlike the picture in Fig. 2, the pi-dependence of the reaction cross section 
displayed for kinematics II in Fig. 7 covers the QFP region up to its boundary, 
i.e., the maximum p\ value is allowed under the condition p r = 0. There has been 
the expectation that the FSI influence is negligible in the QFS region. However, as 
is evident from Fig. 7, the differential cross section is ~ 2 times smaller relative 
to the PWIA predictions in the QFS region due to the FSI effects. It is the 
consequence of the vicinity of the two peaks (QFS and FSI). It is interesting to 
note that the PWIA and FSI predictions for the polarization observables (Figs. 8- 
11) in the QFS region practically coincide while these observables change their 
behaviors very quickly when going to the FSI region. 

6 Conclusion 

In the present paper we have explored a reasonable approximation to treat the 
deuteron-proton breakup reaction in the GeV-region. The employed approach 
incorporates the leading terms in the multiple scattering expansion of the reac- 
tion transition operator in the N — N t-matrix. Our consideration is free from 
shortcomings made in earlier theoretical studies of this process. The differential 
cross sections and the first-order polarization observables have been calculated 
under the two kinematic conditions (one of them corresponding to the Gatchina 
experiment). These calculations have shown that the FSI contributions are very 
important to describe deuteron breakup by the proton in the GeV-region. We 
have seen that the FSI effects, dependent on kinematic conditions, appear in 
different ways. They can slightly change the corresponding cross section with 
unpolarized observables and vice versa. Therefore, studying the deuteron struc- 
ture (especially, the non-nucleonic degrees of freedom in the deuteron) by high- 
energy breakup reactions, one should make every effort to include FSI correc- 
tions. 
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Appendix 

The NN amplitudes Fj in the laboratory frame are related to A-F, the ones in the center-of-mass 
system by 
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. TT (A B\ ^/ U 2 
F 3 = -tU [ — - — I + C ( 1 + 



/C 2 V K X K: 



* = A i& 2 +B + tCU (^-lk 



Fe = D^?z) + £(q*x)}{{fz)--jL(q*x)}+F z) + ±L(Q*x)} {{& z) - 
Fr = D^qx) - ^-(9**)} {(«**) - ^(^as)} +f{(Q*x) - ^(Q**)} |(Q*«) - |^(<2**)} 

|p x pi I _ ppisinOr 



Here we have introduced 



'\2 



(p + p) , s' = (pi +p' 2 ) 
E p +m i _ Epi + E p > 2 



wo = p — Wo = 



75 u ~ 77 



*i = \(u'o + l)(2m + 77) + i(£ pi - £ p , ) 
K 2 = \{u' + l)(2m + 77) - \{E P1 - E p 0. 
The vectors q* and Q* defined in the c.m.s. can be written in the lab. system as 

(q*x) — —(Q*x) = —pisinOi 

J m + ^Ts/2 E P1 + 77/2 \ 



(q*z) = —picos6i+p 



1^5(1+ wo) + 77(1 + U ' ) J 



v ' ^\75(1 + wo) 77(1 + U ' ) 



and their moduli are expressed through the invariants s and t via 



h^-V^f-t, Q* 2 = i(7s + 7s 7 ) 2 +t-4m 2 . 



*2 
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Figures. 

Fig.l Kinematic variables for the pd — > ppn reaction in the deuteron rest frame. 

Fig. 2 The breakup differential cross section at 8i = 15°, 82 = 78° [17) plotted versus the 
fast proton momentum p\ in the deuteron rest frame. The dashed line corresponds to the 
PWIA; the dotted, dash-dotted and solid curves represent the full calculation with J ma x =2, 
3, and 4, respectively. 

Fig. 3 The pi-dependence of the vector analyzing power of the initial proton Co, y , 0,0,0 at 
0i = 15°, 6 2 = 78°. Notations as in Fig.2. 

Fig. 4 The pi-dependence of the deuteron vector analyzing power C V: 0,0, 0,0 at 81 — 15°, 
82 = 78°. Notations as in Fig.2. 

Fig. 5 The pi-dependence of the deuteron tensor analyzing power C,,,, 0,0, 0,0 at 8i = 15°, 
82 = 78°. Notations as in Fig.2. 

Fig. 6 The pi-dependence of the deuteron tensor analyzing power C zz ,o,o,o,o at #1 = 15°, 
82 = 78°. Notations as in Fig.2. 

Fig. 7 The breakup differential cross section at 81 = 6.77°, 82 = 83° plotted versus the fast 
proton momentum pi in the deuteron rest frame. The dashed line corresponds to the PWIA; 
the dash-dotted and solid curves represent the full calculation with Jmax =1 and 2, respectively. 

Fig.8 The same as in Fig. 3 at 81 = 6.77°, 82 = 83°. Notations as in Fig. 7. 

Fig.9 The same as in Fig. 4 at Ox = 6.77°, 82 = 83°. Notations as in Fig. 7. 

Fig.10 The same as in Fig. 5 at 81 = 6.77°, 82 = 83°. Notations as in Fig. 7. 

Fig.ll The same as in Fig. 6 at #1 = 6.77°, 82 = 83°. Notations as in Fig. 7. 
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